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Abstract 

The fundamental structure of the full set of solutions of the BCS 3 P2 pair- 
ing problem in neutron matter is established. The relations between different 
spin-angle components in these solutions are shown to be practically indepen- 
dent of density, temperature, and the specific form of the pairing interaction. 
The spectrum of pairing energies is found to be highly degenerate. 
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Since the discovery of superfluidity in liquid 3 He [[|], great advances have been made 
in understanding the properties of superfluid systems with triplet pairing. In addition to 
the well-studied case of liquid 3 He below 2.6 niK [0-|J, triplet pairing is expected to occur 
in neutron matter in the quantum fluid interior of a neutron star [PHTT||. Neutrino cooling 
processes are strongly affected by 3 P2 pairing in this region |I2"|, as is the vortex structure 



of the star and the coupling between core and crust |l3|Jl^| . A number of common features 



of superfluid, thermodynamic, and magnetic properties of different pair-condensed systems 
have been revealed by analyses based on symmetry principles @-f|, and further analytical 
insights have been gained near the critical temperature T c by application of the Ginzburg- 
Landau approach 0. However, new universalities of triplet pairing may be uncovered by a 
direct attack on the BCS gap equation, as we shall now demonstrate. 

The purpose of this letter is to identify fundamental solutions of the triplet pairing 
problem in neutron matter and elucidate their structure and their relationships. If two 
identical spin-| fermions are paired with a nonzero total momentum J = L + S, the ordinary 
S- wave gap equation is converted into a system of coupled integral equations. In the standard 
notation |S|||, we have the expansion A(p) = J2 ^i!j(p)^i!j( n ) °f the S a P i n the spin-angle 
matrices G^(n; Sl ,s 2 ) = E C C (^ Ms Ml Y LMl ( n) , with (Aft(p))* = (-1) J - M A L f(p) 
assuming time-reversal invariance. The particle-particle interaction has the corresponding 
expansion V(p,p') = E{p\Vi;j J \p)Gfj(n)G^*j(n'), with \L - L'\ < 2 in the case of tensor 
forces. The generalized BCS system then reads [[J 

AiS(p)= E (-^V^^I^^^^O^^tanh^dr', (I) 



VL x J\Mx 



where A = L - L x + 1, dr = p 2 dpdn = dr dn, and S^j^n) = Tr [GfJinjGf^n 
accounts for the summation over spin variables. The energy denominator E(p) = [£ 2 (p) + 
D 2 (p)} 1 / 2 involves the single-particle excitation energy of the normal system and a gap 
function whose square is constructed as 

d\v) = \ E ±%m%jMs%£jilp) ■ (2) 

Z LJML x J\Mx 



The pairing parameter A measuring the gap value at the Fermi surface is given by A 2 = 
/D 2 (pjn)dn/47r. 

Due to the nonlinearity of the gap equation (1), one must in general deal with off-diagonal 
effects in both the total (J, J\) and the orbital (L,L',Li) angular momentum quantum 
numbers. However, for the present we follow the usual practice dating back to Ref. and 
suppress these effects, thus allowing for superposition of spin-angle components only in the 
magnetic quantum number M. The analysis is greatly facilitated by a generalization of 
the separation method developed for S-wave pairing in Ref. fl5f . Thus, defining 4>lj(p) = 



(p\Vlj\Pf)/ v f and Vp = {pf\V^j \pf) , we employ the decomposition 

WtiW) = v f <PlAp)^lAp') + W LJ (p,p') (3) 

of the relevant pairing matrix into a separable portion and a remainder Wlj{p,p') that 
vanishes identically when either argument is on the Fermi surface. Integrals containing Wlj 
as a factor are guaranteed to receive their overwhelming contributions some distance from the 
Fermi surface. In such integrals, the replacements E(p, n) — ► |£(p)| and tanh(£ l /2T) — > 1 are 
justified to high accuracy, errors in neutron matter being of relative order D 2 (p)/e 2 F ~ 1CT 6 , 
where ep is the Fermi energy. 

Substituting (3) into (1) and invoking the orthogonality relation / SffJ^j(n)dn = 5mm 1 , 
the gap equations are recast as 

Bfj = - E / ^(p)5^}(n)^|M tanh ^dr . (5) 

The quantities Bffj are merely numerical factors. Consequently, the p dependence of all gap 
components is seen to be identical. Specifically, we may write A^j(p) = D^jXlj(p), where 
the shape factor Xlj(p) obeys an integral equation 

xlAp) + J w ^M)^y\ K = 0ljOo) (6) 
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of the same form as in the singlet case fL5| . To determine the amplitude -D^j, we note that 
Xlj{pf) = 4>lj{j>f) = 1 since W LJ (p F ,p') = 0. Therefore Afj(p F ) = D^j, and Eq. (4) 
implies Dffj = v F Bffj while Eq. (5) gives 

D% = - VF ^J W«(n)^f^ tanh ^dr . (7) 

The system (6)-(7) is more convenient for solution than the original equations (1), since 
the problem has been divided into (i) evaluation of the M-independent shape factor Xlj(p) 
from the nonsingular linear integral equation (6), and (ii) determination of the structure 
coefficients Dffj from the nonlinear equation (7), where the log-singularity has been isolated. 

Henceforth we specialize to the case L — S — 1, J — 2, this being the most favored 
uncoupled channel for pairing in neutron matter at densities prevailing in the quantum 
fluid interior of a neutron star (k F = pp/h ~ 2 fm _1 ), where the 1 So gap has already 
closed P,|15[- The arguments are simplified if we adopt D® 2 = <5 as a scale factor, write 
Di2^° = (Am + iKAi)fi/V6) an d introduce a "structure function" 

d 2 (n) = l67rD 2 (p)/xUp) = 5 2 [(1 + A 2 ) 2 + «? + «£ + (\\ - 4A 2 - k\)x 2 
- 2(Ai + A X A 2 + k x k 2 )xz + (3 + \\ - \j - 2\ 2 )z 2 

+ 2(2k 2 - KiXx)xy + 2(«i + X 1 k 2 - \ 2 K 1 )yz] , (8) 

where x = sin # cosy?, y = sin^sinyj, and z = cos 8. After separation of the real and 
imaginary parts of the -Df 2 ^° in Eq. (7), we arrive at the set of equations 

A 2 = -v F [X 2 (Jo + Js) ~ Ai^i - KiJ 2 - J3] , 
n 2 = -v f [k 2 {Jq + J 5 ) - K\J\ + Ai J 2 + J4] > 
Ai = —v F [\i Jq - (A 2 + 1) Ji + k 2 J 2 - k 1 J 4 /2] , 

Ki = —Vp^Jj - K 2 Ji - (A 2 - 1) J 2 - Ai J4/2] , 

1 = -v F [-(X 1 Ji - ki J 2 + A 2 J 3 - «2-/ 4 )/3 + J 5 ] , (9) 

which, in angular content, is consistent with the corresponding set in Ref. ||. The integrals 
J k are given by J 6 = ( J + 4J 5 + 2J 3 )/4, J 7 = ( J + 4J 5 - 2 J 3 )/4, and, for A; = 1, • • • , 5, by 
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J k = [ue^ f^fU ^^dr, (10) 

J J K 2E(p) 2T K J 

with /o = 1 - 3z 2 , h = 3xz/2, f 2 = 3yz/2, f 3 = 3(2x 2 + z 2 - l)/2, / 4 = 3xy, and 
f 5 = (l + 3z 2 )/2. 

The system (9) has three one-component solutions |8||| with \M\ = 0, 1, and 2. We 
preface our analytic exploration of multicomponent solutions with the following observation. 
Substitution of dd 2 (9 , ip) / dip for f k in definition (10) must yield zero upon integration over 
ip. This identity implies a relation 
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£ Cfe j fc = o ( n ) 

fc=l 

between the integrals, with c\ = K\ + A1K2 — A2K1, C2 = Ai + A1A2 + ^1^2, C3 = 2k 2 — KiAj, 
and c 4 = 2A 2 — (A^ — nl)/2. Now observe that if the first equation of (9) is multiplied by 2k 2 , 
the second by 2A2, the third by K\, and the fourth by — Ai, and the results of the last three 
operations are subtracted from that of the first, relation (11) is reproduced. Thus only four 
of the five equations in (9) are truly independent and hence any one of the parameters Ai, 
A2, Ki, K 2 can be chosen arbitrarily. We take K\ = 0. With this choice, solutions of (9) are 
necessarily even functions of Ai, so attention may be focused on the sector Ai > 0. 

The search for multicomponent solutions begins with the restricted case k 2 = 0, for 
which d 2 (n) is independent of y. In Eq. (10), the integration of fk(0,<p) over y is then 
carried out using the formula sm 9 d9 dip = 25(x 2 + y 2 + z 2 — l)dx dy dz. For k = 1,3, this 
yields 2fk(x, z)(l — x 2 — z 2 )~ x l 2 since fi and ^3 are independent of y, while both J 2 and 
J4 = vanish because f 2 and f& are odd in y and the y integral has symmetric limits. As a 
result, there remain only three independent equations, 

\ 2 = -v F [\ 3 (J + J B )-\iJi-J 3 ], (12) 
Ax = -uf[Ai(J /4 + J 5 ) - (A 2 + 1)J X + A1J3/2] , (13) 
1 = -w F [-(AiJi + A 2 J 3 )/3 + J 5 ] . (14) 

We first identify and verify a particular solution with Ai = 0, A 2 = 3, for which the structure 
function (8) becomes d 2 (x, z) = A5 2 [A — 3(x 2 + z 2 )}. The symmetry of this function with 
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respect to x and z implies the relation 3Jo + 2J 3 = 0, since the combination 3/o + 2/ 3 = 
6 (a; 2 — z 2 ) changes sign on interchange of x and z whereas d 2 and other factors within the 
integrand of (10) are left unchanged. Further, Ji(Ai = 0, A2 = 3) = since fi = xz is an 
odd function of x. Under these conditions, Eq. (13) is satisfied identically, while Eqs. (12) 
and (14) coincide and the resulting equation, 1 = —vp(J 5 — J 3 ), determines 8. All other 
solutions of the set (9) are more degenerate. To illustrate this important feature, let us put 
A 2 = —1. Then, at any Ai the structure function (8) is seen to take the factorized form 
d 2 {n; Ai, A 2 = -1) = 5 2 (A 2 + 4){x 2 + z 2 ) = 24vrA 2 (x 2 + z 2 ). The symmetry of d 2 in x and 
z again implies the relation 3Jo + 2J 3 = 0, while integration of f\ = xz over x gives and 
therefore J\ = 0. It follows that Eqs. (12)— (14) again coincide but now provide an equation 
1 = — vjr[J 3 /3 + J 5 ] that determines A 2 rather than A x or 5 individually. Here we have 
a striking example of the universal structure of solutions of the 3 P 2 pairing problem, also 
manifested in the remaining solutions of the system (9). 

These further solutions are found by implementing a rotation R — (x — tcosa + 
usma, z = t sin a — mcosq;). Expressing d 2 in terms of t and u and setting tana = 7, 
one easily finds conditions 

(A 2 -4A 2 ) 7 2 + 2A 1 (l + A 2 ) 7 + A 2 -A 2 -2A 2 + 3 = 0, A l7 2 - (A 2 - 3)7 - A x = (15) 

under which d 2 becomes a function of t only. The choice 7 (Ai, A 2 ) = 7o(Ai, A 2 ) = 
+ A 2 )/(4A 2 — A 2 ) meets both conditions provided 

(A 2 -2A 2 + 2)(A 2 -2A^-6A 2 ) = 0. (16) 

Equation (16) embodies three branches of A 2 versus Ai, which start as parabolas from Ai = 
and A 2 = 1, 0, and —3. The structure function has t dependence d 2 (t) oc 1 — t 2 when the 
first factor of (16) vanishes and d 2 (t) oc 1 + 3t 2 when the second is zero. Calculating the 
integrals J\ and J 3 by rotation of the x, z plane under R, we are led to the relations 

(A 2 -4A 2 )J 1 + A 1 (A 2 + 1)J 3 = 0, 
3Ai(l + A 2 ) J - 2(A 2 - 2A 2 + 6) h = . (17) 
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The first relation (for example) is verified as follows, noting that the integrand on its l.h.s. 
is proportional to [(A 2 - 4A 2 )/i + Ai(A 2 + l)/ 3 ]/(l - t 2 - u 2 ) 1/2 . Substituting fi(t,u) and 
fs(t,u) and integrating over u, which can be done freely for any shape of d 2 (t), we obtain a 
result that is proportional to (Af — 4A 2 )7 + Ai(A 2 + 1) and therefore vanishes when 7o(Ai, A 2 ) 
is substituted. What is remarkable is that Eqs. (12)-(14) coincide when relations (17) are 
inserted, and once again these equations determine only A 2 . Thus, construction of the 
rotation R "kills two birds with one stone": the condition (16) required to transform d 2 
into one-dimensional form also specifies another set of solutions of our system. Within the 
constraint K\ = 0, these solutions possess a line degeneracy as opposed to the point (or 
nondegenerate) character of the solution (Ai = 0, A 2 = 3). In addition to the free choice 
made for k±, one of the coefficients A« can be assigned arbitrarily. 

We now allow k 2 to have a nonzero value, thus bringing in the second of Eqs. (9). At 
Ai = 0, this equation becomes identical with the first of the set, as is seen with the aid of 
Eq. (11). The particular solution (Ai = 0, A 2 = 3) is then replaced by one with Ai = and 
(A 2 . + ft 2 ) 1 / 2 = 3, but the pairing energy remains unaltered, the relevant quantities being 
independent of the phase of the coefficient -D 2 2 (Ai = 0). To find the other multicomponent 
solutions in the general case with n 2 ^ and Ai ^ 0, we may extend our previous tactic and 
apply a rotation in t/iree-dimensional space so as to eliminate four terms in expression (8) 
and cast d 2 into a one-dimensional form. The three Euler angles are thereby fixed, implying 
the single relation 

4 = (l + A 2 )(A 2 /2- A 2 + 1) (18) 

between A 2 , re 2 , and Ai, which can be shown to satisfy all of Eqs. (9). This relation defines 
two branches « 2 (Ai, A 2 ). Starting at the plane k 2 = 0, one branch grows out of the solution 
A 2 = —1 while the other grows out of the parabola A 2 /2 = A 2 — 1 contained in Eq. (16). 
(Accordingly, A 2 = — 1 and this parabola cannot be counted as independent solutions.) The 
two surfaces defined by (18) complete the set of states of the 3 P 2 problem. 

The solutions we have identified divide into two groups, the states within a group being 
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essentially degenerate in energy. This behavior is consistent with the numerical calculations 
reviewed in Ref. ||. The group with lowest energy, having structure function d 2 (t) oc 1 + 3t 2 , 
contains only nodeless states and consists of (i) the particular state (Ai = 0, \J X 2 . + k 2 = 3) 
and (ii) the states belonging to the branches of Eq. (16) starting at the points (Ai, A2) = (0, 0) 
and (0, —3). The upper group contains the remaining states, having d 2 (t) oc 1 — t 2 and one 
node. The splitting between the two groups can be calculated (for example) as the splitting 
between the states A2 = —1 and (Ai = 0, A2 = 3), henceforth labeled u and 1 respectively. 
Evaluating (14) for both states, one finds the relation J3 /3 + J3 + Jg — jj = between 
integrals of the form (10). In explicating this relation, we exploit the fact that the dominant 
contributions to the integrals J3 , and — J^ n come from the range of p values 
adjacent to the Fermi surface, where Xvzip) an d 0ia(p) are effectively unity. One readily 
arrives at the analytical result 

"'i^ '^^ 1 " 3 "- - 028 (19) 

for the splitting of upper and lower states, in close agreement with Ref. ||. Similar results 
are also available at finite temperature T. 

The conclusions that follow from these exercises are that if the mixing of different L, J 
channels is neglected, (i) the 3 P2 gap spectrum is nearly degenerate and (ii) its structure, 
in terms of energy splittings between the different states, is a universal function of T/T c , 
independent of any other input parameters including the density. In particular, the concrete 
form of the particle-particle interaction V was not used anywhere, so the structure and 
relations we have established retain their validity even when fluctuation and polarization 
corrections to the bare V are taken into account. 

Finally, we return to the issue of nondiagonal contributions to the system (1) of gap 
equations, which arise principally from the 3 P2- 3 F 2 coupling, and outline a perturbative 
evaluation [[| of their effects. The r.h.s. of each equation of the set (9) is now perturbed by a 
small "nondiagonal" contribution. In the presence of these additional terms, the degeneracies 
found above are removed. The pairing energy no longer depends on d 2 (n) alone, and the 
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parameters Aj, Ki are fully determined. Consider, for example, the alteration of the parabolic 
branch contained in the relation (16), which may be measured by a new variable C, = 
XI — 2X2 + 2. Both ( and the change 77 = 5(() — 5(( = 0) of the scale factor 5 are expected to 
be small; therefore in performing Taylor expansions of the Jk((,rj,\i) we need only retain 
terms linear in ( or 77. The original set of three equations (12)-(14) is replaced by three new 
ones, each of which takes the schematic linear form (A(Xi) +r]B(Xi) = -P(Ai) with different 
choices of the functions A, B, and P, all referred to ( = 77 = 0. The small quantities ( and 
77 may be obtained from any pair of the equations, as functions of Ai. Substituting these 
functions into the remaining equation, we arrive at a closed form that determines the value 
of Ai, which was hitherto arbitrary. Estimation and analysis of available numerical results 
IHJ indicate that the nondiagonal corrections to the universal relations we have derived 
for triplet pairing in neutron matter are small, maximally of order several percent of the 
splitting given by Eq. (19). 



In summary, straightforward arguments based on a new separation method |TjJ for treat- 
ing BCS-type gap equations have revealed the structure and energetics of the full set of 
solutions of the pairing problem in the uncoupled 3 P2 channel. In contrast to the Ginzburg- 
Landau scheme employed by Mermin ||, the present approach is applicable at any tem- 
perature T. The analysis shows that the structure of the solutions is in fact universal - 
independent of the temperature, the density, and the specific parameters of the interparticle 
potential, which affect only an overall scale factor in the pairing energies. The line of analysis 
we have followed transcends the problem considered here. An obvious future objective is to 
characterize the solutions occurring in the problem of superfluid 3 He. The structure function 
d 2 (x,y,z) is again bilinear in its variables, but the states for L = S = 1 with J = 0, 1,2 
contribute on an equal footing and the number of equations in the system analogous to 
(9) rises from five to nine. The treatment introduced herein could also be relevant to the 
description of superdeformed bands in atomic nuclei, if triplet P-wave pairing is responsible 
for this phenomenon as suggested in Ref. [|T5]]. 
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